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Abstract 

In this paper, we are concerned with the structural stability of some steady subsonic 
solutions for Euler-Poisson system. A steady subsonic solution with subsonic background 
charge is proven to be structurally stable with respect to small perturbations of the 
background charge, the incoming flow angles, the normal electric field and the Bernoulli's 
function at the inlet and the end pressure at the exit, provided the background solution 
has a low Mach number and a small electric field. Following the idea developed in [19] . 
we give a new formulation for Euler-Poisson equations, which employ the Bernoulli's law 
to reduce the dimension of the velocity field. The new ingredient in our mathematical 
analysis is the solvability of a new second order elliptic system supplemented with oblique 
derivative conditions at the inlet and Dirichlet boundary conditions at the exit of the 
nozzle. 
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1 Introduction 



The Euler-Poisson (or hydrodynamical) model for a unipolar semiconductor in the isen- 
tropic steady state case reads 



(P u i)x! + {pUlU 2 )x 2 + (pUlUs) X3 +p Xl = ptp xi . 



{pUiU 2 ) Xl + {pu\) X2 + (P u 2U3)x s +Px 2 = P<P 
(puiu 3 ) Xl + (pu 2 u 3)x2 X3 — P*P 



X2 - 



X-.i i 



A(p = p — b{x). 



(1.1) 
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where p, u, ip denotes the electron density, the average electron velocity and the electrostatic 
potential, b(x) be the prescribed ion background density (doping profile). The electric field 
E is given by E = V</?. The energy equation of the hydrodynamic model is replaced by the 
pressure-density relation p = p{p). One may refer to |14l [T6] for more details on physical 
background. 

In this paper, we are concerned with the structural stability of some steady subsonic 
solutions for Euler-Poisson system. Indeed, we will construct a subsonic solution to (II. ip in a 
rectangular cylinder by imposing suitable boundary conditions, which is also required to be 
close to some special subsonic solutions. The boundary conditions we will impose have their 
origins on our previous works in [9l IT71 [T8| [T9] and another important work [12] on transonic 
shock solutions for (jl.ip . As we have done in [19] , we will prescribe the incoming flow angles 

fa = —,@3 = — and Bernoulli's function B = -fit? + Mo + u\) + hip) — ip at the inlet, the 
u\ u\ 2 

natural slip boundary conditions on the nozzle walls and the end pressure at the exit. For 
the electric field E = V<^, we prescribe the normal component of E at the inlet and nozzle 
walls, and ip at the exit of nozzle. 

The one-dimensional steady state isentropic flow hydro dynamical model was analyzed in 
[7J and the three-dimensional irrotational case was disucssed in [8]. In both papers, existence 
and uniqueness results for small data generating subsonic flow were proven. In |15j . the 
author obtained the 2-D existence of smooth solution under some smallness assumptions on 
the prescribed outflow current and the gradient of the velocity relaxation time. Also, the 
author used a smallness assumption on the physical parameter multiplying the drift-term 
in the velocity equations instead in the irrotational 3-D case and obtained similar results. 
For the full hydrodynamic model,Yeh [20] showed the existence of a unique strong solution in 
several space dimensions if the flow is subsonic, the ambient temperature is large enough, and 
the vorticity on the inflow boundary and the variation of the electron density on the boundary 
are sufficiently small. Zhu and Hattori [21] proved the existence of classical subsonic solutions 
in one space dimension for the whole space problem under the additional assumption that the 
doping profile be close to a constant. See also [4 J and the reference therein for more details. 

Recently, there are some important progress on transonic shock solutions for 1-D Euler- 
Poisson systems \\.2\ I13j. Gamba [IOj showed existence of steady-state solutions in the tran- 
sonic case by means of the vanishing viscosity method. However, the solutions as the limit of 
vanishing viscosity may contain boundary layers and more than one transonic shock. In [12] , 
Luo and Xin gave a thorough study of the structure of the solutions to boundary value prob- 
lems for 1-D Euler-Poisson system for different situations when the density of fixed, positively 
charged background ions is in supersonic and subsonic regimes. The existence, non-existence, 
uniqueness and non-uniqueness of solutions with transonic shocks were obtained according 
to the different cases of boundary data and physical interval length. The solutions they 
constructed contained exactly one transonic shock in the interval [0, L]. Moreover, they can 
determined the shock location by the boundary data and L. In [13] . the authors investigated 
structural and dynamical stabilities of steady transonic shock solutions for one-dimensional 
Euler-Poisson systems. It was shown that a steady transonic shock solution with a supersonic 
background charge was structurally stable with respect to small perturbations of the back- 
ground charge, provided that the electric field is positive at the shock location. Furthermore, 
any steady transonic shock solution with a supersonic background charge was proved to be 
dynamically and exponentially stable with respect to small perturbations of the initial data, 
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provided the electric field is not too negative at the shock location. The remaining natural 
question is the structural stability of these transonic shock solutions in a multi-dimensional 
domain supplemented with suitable boundary conditions, a considerably more difficult mat- 
ter. We start to investigate the structural stability of some special subsonic solution to take 
a close look at the physically acceptable boundary conditions we should prescribe and hope 
this may serve as a building block toward that challenging goal. 

In 01171118], we have characterized a set of physically acceptable boundary conditions that 
ensure the existence and uniqueness of a subsonic irrotational flow in a finitely long flat nozzle 
by prescribing the incoming flow angle and the Bernoulli's function at the inlet and the end 
pressure at the exit. In [19], we have developed a new formulation for the three dimensional 
Euler equations. The key idea in our formulation is to use the Bernoulli's law to reduce the 
dimension of the velocity field by defining new variables (l,fa = — ~,fa = ~~~) an d replacing 

u% by the Bernoulli's function B. We find a conserved quantity for flows with a constant 
Bernoulli's function, which behaves like the scaled vorticity in the 2-D case. Moreover, a 
system of new conservation laws can be derived, which is new even in the two dimensional 
case. Following the ideas developed in [19], we reformulate the Euler- Poisson equations in 

terms of (s = lnp, fa, /3s, B,(p) and replace u\ by B through u\ = ~!L aT — ' Here 

one should note that due to the Bernoulli's law (|2.2p . the Bernoulli's function B will possess 
the same regularity as the boundary data at the inlet, which is quite different from the 
velocity field. In this way, we can explore the role of the Bernoulli's law in greater depth 
and hope that may simplify the Euler equations a little bit. We can also find a new quantity 
W = 82^3 — 83 fa + fadifa — fadifa, which is conserved along the particle path at least for 
flows with a constant Bernoulli's function (see (|2.6p ). In subsonic region, we find that the 
density and the potential function are coupled together to satisfy an second order elliptic 
system, which will possess good regularity. Roughly speaking, the equation (|2.6p shows that 
W possess one order lower regularity than those of the density and the potential function. 
This, together with the first three equations in (12.3P form an divergent-curl system for (fa, fa), 
which indicates the velocity field (fa, fa) should possess the same regularity as those of the 
density and the potential function. 

We make some comments on our proof. The new ingredient in our mathematical analysis 
is the solvability of some second order elliptic system supplemented with oblique derivative 
conditions at the inlet and Dirichlet boundary conditions at the exit of the nozzle. Indeed, 
both p and ip satisfy some second order elliptic equations and they are coupled together to 
form a second order elliptic system. Besides this, the boundary conditions for p and (p are 
also coupled together at the inlet of the nozzle. This elliptic system is neither weakly-coupled 
nor cooperative in the sense of p]. To guarantee the uniqueness, some smallness assumptions 
are prescribed on the background solution. That means our background solution should have 
a low Mach number and a small electric field. 

This paper proceeds as follows. The second section contains a new formulation for Euler- 
Poisson system and a statement of our background solution. The last section focus on the 
structural stability of our background solution with respect to small perturbations of the 
background charge, the incoming flow angles and the end pressure, provided the background 
solution has a low Mach number and a small electric field. 
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2 Preliminary 



2.1 A new formulation for Euler-Poisson equations 

Following [19], we develop a new formulation for Euler-Poisson equations, which employ 
the Bernoulli's law to reduce the dimension of the velocity field. By employing the following 
identity in vector calculus 

u • Vu = V(-|u| 2 ) — u x curlu, 
and the momentum equation in (jl.ip . we have 

V(^|u| 2 + h(p) - ip) - u x curlu = 0. (2.1) 
Then the Bernoulli's law holds: 

u • VB = 0. (2.2) 

Here the Bernoulli's function is defined to be B = — (uf + u\ + u|) — tp + h(p), where h'(p) = 

c 2 (p) = p'(p) 
P P 

U 2 U 3 

Define the following three new variables: f3 2 = — , /3 3 = — ,s = In p. Then u\ has a 

Ml U\ 

. 2 2{B + <p-h{p)) 
simple expression: u-i = -s -s — . 

l + /3 2 2 + /?| 

— ^ * 1 

Multiplying the first equation in (11.11) by ^5— , dividing the i-th equation in (jl.ip by 

puf 

pu\ and adding them together for i = 2, 3, 4, we obtain the following new system: 

c 2 (p) 1 

#is + /3 2 92« + /33<9 3 s + d 2 /3 2 + d 3 f3 3 = — =di<p, 

uf uf 

difo + fod 2 p 2 + /W2 - ^/3 2 5is + ^-d 2 s + - = 0, 

uf uf uf uf 

dxfc + /w 3 + /Ws - ^/? 3 a lS + c ^d 3 s + %d x y - = o, ( 2 - 3 ) 

uf uf uf uf 

diB + p 2 d 2 B + (3 3 d 3 B = 0. 
= p — b(x). 

Then it is easy to show s and (p will satisfy the following second order elliptic system in 
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the subsonic region, i.e. |u| < c(p), 

dx ( - s - fod 2 s - Mss) +d 2 ( - fodis + - $)dzs - f3 2 (3 3 d 3 s 



-ih ! - hih -s - ft/W + - Z^f)^ - f (4 " ^ - - f3 3 8 3 s 



+ ( (9 2 /3 2 ) 2 + (<9 3 /3 3 ) 2 + 2«9 2 /W 2 J - 4( e " " "V^vi = 0, 

= e s - b(x). 

(2.4) 

An important observation made in [TU] is the quantity W = d 2 f5 3 — d 3 j3 2 + f3 3 di/3 2 — f3 2 d\f3 3 
is conserved along the particle path for flow with a constant Bernoulli's function. To simplify 

the notation, we set G = l+/3f +/3f and D = di+fi 2 d 2 +p 3 d 3 . Plugging u\ = 2 ^ B + J 2 ~ h ^ 

1 + p 2 + P3 

into the second and third equation in (|2.3p . then we obtain 



DA - 2( B + r -M P )) ° (ftaiS " * S) + 2(B + ff - M p)) (&8iy - **> - (2 . 5) 

Apply j3 3 d\ — <9 3 and —fi 2 d\ + ^2 to the above two equations respectively and add them 
together, one can show that W satisfies the following equation. 

D © + Mwrh m? ( (1> *■ ft) ■ I(£ r v " " Vv) x VB| ) = a <2 ' 6) 



The verification of (|2.6p is similar to the calculation in the appendix in [19]. For com- 
pleteness, we also give a detailed calculation for (]2.6p in the appendix. 
Suppose that B = const, then (j2.6|) reduces to 

D(^)=0. (2.7) 

W 

This implies — is conserved along the particle path. Indeed, since B = const, by (12. lh . 

the vorticity field is parallel to the velocity field. So we may assume that there exists a real 
function p{x) such that curlu = a(x)u. By simple calculations, we have 

W curlu • u ^(x)|u| 2 fj,(x) (2 8) 

~pG ~ pu\G ~ pu\G ~ ~p~' 

Suppose that B = const, then (j2.6|) reduces to 

D (^)=0. (2.9) 

Historically, the stationary solution of Euler equations with the vorticity field paralleling to 
the velocity field was called Beltrami flow and have been investigated for over a century. One 
may refer to Arnol'd [2], Constantin and Majda [3] for more details. Here we emphasize that 
our calculations works for any general Euler-Poisson flows. 
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Remark 2.1. As we have discussed in the introduction, W may help to raise the regularity of 
the velocity field fii and P3. In subsonic region, the equation {2.1$ shows that the density and 
the potential function are coupled together to satisfy an second order elliptic system, which 
will possess good regularity. Roughly speaking, the equation \2. 6\) shows that W possess one 
order lower regularity than those of the density and the potential function. This, together with 
the first three equations in \2. 3\) form an divergent-curl system for (,$2, /3a), which indicates 
the velocity field (/32j/?3) should possess the same regularity as those of the density and the 
potential function. 

Remark 2.2. It is easy to see that G satisfies the following Riccati-type equation: 

DG " 1 U ^ G2 di(%) - <P) - ,„ , 1 h( ss GT>(h(p) -<p) = 0. (2.10) 

(B + (p-h{p)) {B + (p-h{p)) 

One may expect some blow-up results. However, the ambiguous sign of d\(h(p) — ip) in the 
coefficient of G 2 in t2.10\) makes the whole argument nontrivial. One should note that G 
blows up means that the fluid will turn around in the flow region. 

Remark 2.3. In U9\j , we can also derive a system of new conservation laws for compressible 
Euler equations. However, we can not obtain similar conservation laws for the Euler-P oisson 
equations due to the effect of the electrostatic potential. 

Remark 2.4. p and ip are coupled together to form a second order elliptic system. This elliptic 
system is neither weakly- coupled nor cooperative in the sense of JiJ/. Hence the uniqueness 
for {2.1$ with suitable boundary conditions will be a main obstacle in our analysis. 



2.2 Background solutions 



(2.11) 



In |12j . the author consider the initial value problem for the following 1-D Euler-Poisson 
equations 

(pu) x = 0, 
(jp(p) + pu 2 ) x = pE, 
E x = p- b . 

(p,u,E){0) = (p I ,u I ,E I ). 

Here 60 is a positive constant. 

Assume uj > 0. By the first equation in (|2.1ip . we have pu(x) = J = pjuj and the 
velocity is given by 

u = J I p. (2.12) 
Thus the boundary value problem (|2.11|) reduces to 



j2 

(P(.P) + — )* = pE, 
P 

E x = p- b , 



(2.13) 



(p,E)(0) = (p I ,E I ). 

Use the terminology from gas dynamics to call c = y/p'(p) the sound speed. There is a 
unique solution p = p s for the equation 

P(p) = J 2 /P\ (2-14) 
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which is the sonic state (recall that J = pu). Later on, the flow is called supersonic if 

p'(p) < J 2 /p 2 , i.e. p < p s . (2.15) 

Similarly, if 

p\p) > J 2 /p 2 , i.e. p > p s , (2.16) 

then the flow is said to be subsonic. 

The solution of (|2. 13[) was analyzed in (/?, E) — phase plane. Any trajectory in (p, E)— plane 
satisfies the following equation, 

d (lE 2 -H(p)) = 0, where H'(p) = ^V(p - ^)). (2.17) 

L p p 

The trajectory passing through the point (pi,Ei) with pj > is given by 

\E 2 - J*H'{ S )ds= l -E 2 . (2.18) 
We only consider the case b$ > p s , i.e. &o is i n subsonic region. 

Definition 2.5. The critical trajectory (for the case bo > p s ) is the trajectory passing through 
the point (feojO) with the equation: 



\e 2 - f H'(s)ds = 0. (2.19) 
2 Jbo 



We are concerned with the following two cases, that is case cl) and c4) in the subsection 
2.2 in [12]. Suppose (pi,Ej) is on the critical supersonic trajectory, i.e. 



1 f 

-Ej+ / H'(s)ds = 0. 

2 J pi 

p s < pi < b ,E! > 0. 



In this case, initial value problem (12. 13ft admits a unique subsonic solution (p, E) for 
all x > 0. Moreover, 

p x > 0,E X < 0,x > 0, lim (p,E)(x) = (b ,0). (2.20) 

X— >oc 

pi > b ,Ei < 0. 

In this case, initial value problem f)2. 13j) admits a unique subsonic solution (p, E) for 
all x > 0. Moreover, 

p x > 0, E x > 0, x > 0, lim (p, E){x) = (b , 0). (2.21) 
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Due to our technical reasons, we need to choose some special background solutions. Given 
60 > 0, we can choose J > such that 60 > p s , i.e. 60 is in subsonic region. We take pi be 

1 n~ 

close to bo , so that < uj < — * — — - — (we can do it by choose a small J) . Then we choose 

3 y b Cn e 

1 

< \Ej\ < — g and (pj,Ej) is on the critical supersonic trajectory. By the above 

analysis, we can find a unique subsonic solution (pg, Eo) to f|2. 13|) for all x > 0. Moreover, we 
have lim (p,E)(x) = (60, 0). Restricted to [0, 1], (po,E ) will satisfy the following properties: 

x— >oo 

^ 77 < rain p (x) < max Po(x) < ^6 . 
2 xe[o,i] x-e[o,i] 2 



l\ L 1 < min u (x) < max u (x) < - J \ . (2.22) 
6 y b Cn e xe[o,i] xe[o,i] 3 y b Cn e 

max E (x) < „ . 

xe[o,i] C(b )C% e 

Here C(6q) is a smooth function of bo. The constant Cn e is the least number such that the 
following inequality holds: 

\\U(x)\\ 2 L2{ne) < Cn e \\VU(x)\\l 2[ne) ,y U(x) ei= {U(x) G H\n e ) : U(l,x 2 ,x 3 ) = 0.}. 

For convenience, we introduce the electrostatic potential ipo(x) satisfying ^o( x ) = ^o(^) 
with 93(1) = 0. 

Remark 2.6. The background solution we have chosen has a low Mach number and a small 
electric field in the sense of (E 



3 Structural stability of background solutions 



In this section, we are concerned with the structural stability of our background solutions 
for Euler-Poisson system. Indeed, we will construct a subsonic solution to (jl.ip in a rectan- 
gular cylinder by imposing suitable boundary conditions at the inlet and exit, which is also 
required to be close to our background solutions (po(xi),uo(xi),ipo(xi)). The rectangular 

cylinder will be f2 = [0, 1] x [0, 1] x [0, 1]. We also set B = -ul(xi) + h(po(xi)) — (po(xi) is 

a constant, sq(xi) = lnpo(a;i). 

At the inlet of the nozzle x\ = 0, we impose the flow angles and the Bernoulli's function: 

Pi{0,x 2 ,x 3 ) = ePf n >(x 2 ,x 3 ),i = 2,3, 

(3-1) 

B(0,x 2 ,x 3 ) = B Q + eB in {x 2 ,x 3 ). 
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Here the compatibility conditions should be satisfied: 

di(3?(0,x 3 ) = di/3 l 2 n (l,x 3 ) = 0, 
^ n (x 2 ,0) = ^ n (x a ,l) = 0, j = 0,2, 

d k 2 B m {^x 3 ) = dlB in (l,x 3 ) = d$B in (x 2 ,0) = d$B in (x 2 , 1) = 0, k = 1,3. 
At the exit of the nozzle x\ = 1, we prescribe the end pressure: 



(3.2) 



1 



,"/(so(l)+€S e (x2,X3)) 
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p{l,x 2 ,x 3 ) 

Here we also require that p e satisfies the following compatibility conditions: 

^s e (0,x 3 ) = ^s e (l,x 3 ) = 0, 
d J 3 s e (x 2 ,0) = d J 3 s e (x 2 , 1) = 0, j = 1,3. 
While on the nozzle walls, the usual slip boundary condition is imposed: 

u 2 (xi,0, x 3 ) = u 2 (x 1 , l,x 3 ) = 0, 
u 3 (x 1 ,x 2 ,0) = u 3 (x 1 ,x 2 , 1) = 0. 



(3.3) 



(3.4) 



(3.5) 



For the electric field E = Vys, we impose the normal component of E at the inlet and the 
nozzle walls, while at the exit, we prescribe the Dirichlet boundary condition for tp: 



dM0,x 2 ,x 3 ) = E I + eE in (x 2 ,x 3 ), 
d 2 (p(xi,0,x 3 ) = d 2 ip(xi, l,x 3 ) = 0, 
d 3 ip( Xl ,x 2 ,0) = d 3 ip(xi,x 2 , 1) = 0, 
(p(l,x 2 ,x 3 ) = 0. 

Here E m (x 2 ,x 3 ) should satisfy the following compatibility conditions as B m : 

dlE m (0,x 3 ) = dlE in (l,x 3 )=dlE in (x 2 ,0)=dlE m (x 2 ,l)=0, k = 1,3. 



(3.6) 



(3.7) 



The prescribed ion background density b{x) is taken to be b(x) = bo + eb{x). Here 
b(x) € C 1,a (0) should satisfy the following compatibility conditions: 



d 2 b(xi,0,x 3 ) = d 2 b(x 1 , l,x 3 ) = d 3 b(x 1 ,x 2 ,0) = d 3 b(x 1 ,x 2 , 1) = 0. 



(3.8) 



Mathematically, we are going to prove that ([Lip with boundary conditions (|3. 1|) . (|3.3|) . 
(|3.5p and (|3.6p satisfying compatibility conditions ([3.2p , (|3.4p and (|3.7p , has a unique subsonic 
solution. 
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3.1 Extension to the domain Q e = [0, 1] x T 2 



Suppose the flow (p, Ul ,u 2 ,u 3 ) G C 3 ' a (n) x C 2 ' a (n) 3 we will construct has the following 
properties: 



d J 2 (p,u 1 ,tp)(xi,0,x 3 ) = di(p,u lj ip)(xi,l,x 3 ) = 0, 
d 3 i (p,u 1 ,ip)(x 1 ,x 2 ,0) = dl(p,u 1 ,ip)(x 1 ,X2, 1) = 0, j = 1,3, 
<9 2 u 2 (xi,0,x 3 ) = d 2 u 2 (x!,l,x 3 ) = 0, 
dlu 3 ( Xl ,x 2 ,0) = dlu 3 (x u x 2 ,l) =0, k = 0,2. 



(3.9) 



Then we may extend (p, ui, u 2 , 143) in the following way, to (p, iii,u 2 ,u 3 ) G C ,a ([0, 1] x 
C7 2 ' Q ([0,1] x M 2 ) 3 : 

For (x\,x 2 ,x 3 ) G [0,1] x [—1,1] x [—1,1], we define (p,ui,u 2 ,u 3 ) as follows 



(p,Ul,U 2 ,U 3 ,0)(x) 



(p,ui,u 2 ,u 3 ,ip)(x 1 ,x 2 ,x 3 ), if (x 2 ,x 3 ) G [0,1] x [0,1], 

-n 2 ,n 3 , -x 2 ,£ 3 ), if (x 2 ,x 3 ) G [-1,0] x [0, 1], 

(p,u 1 ,u 2 ,-u 3 ,(p)(x±,x 2 ,-x 3 ), if (x 2 ,x 3 ) G [0,1] x [-1,0], 
( / f3,ui,-n 2 ,-ii 3 ,93)(xi,-X2,-x 3 ), if (x 2 ,x 3 ) G [-1,0] x [-1,0]. 



Then we extend (p, , tp),ui,u 2 ,u 3 periodically to [0,1] x M 2 with period 2. It is easy to 
verify that (p, , $), u 1: u 2 , u 3 ) will belong to C 3 ' a ([0, 1] x M 2 ) 2 x C 2 ' Q ([0, 1] x M 2 ) 3 . More- 
over, (p,(p,ui,u 2 ,u 3 ) will also satisfy Euler-Poisson equations. Due to these reasons, one 
may directly work on the domain [0,1] x T 2 (Here T 2 is a 2-torus), so that the difficulty 
caused by corner singularity and slip boundary conditions will be avoided. We may extend 
{Pi n ,Pi n ,B in ,E in ,s e ) to T 2 , which will still be denoted by ($ n , /3| n , B in , E in , s e ). 

3.2 Main results 



The main result is the following existence and uniqueness theorem. 

Theorem 3.1. Given {^ 2 n ,^ 3 n ,B in ,E in ,s e ) G C 3 ' a (T 2 ) and b{x) G C l > a (£l e ), there exists a 
positive small number eo, which depends on the background subsonic state (po, no, 0, 0, (fo) 
and (/3 2 n , /3 3 n , B m , E m , s e , 6), such that if < e < eo, then there exists a unique smooth 
subsonic flow (ui,u 2 ,u 3 , p, <p) G C 2 ' a (f2 e ) 3 x C 3 ' a (tt e ) 2 to (jl.ip satisfying boundary conditions 
(|3.ip - ()3.5p and ()3.6p . Moreover, the following estimate holds: 



\\(u 1 ,u 2 ,u 3 ) - (■u ,0,0)|| c -2, a( n e ) + \\p - p \\c3, a( n e ) + \\<f ~ yo||c3,<*(n e ) < Ce. (3.10) 
Here C is a constant depending on (p Q , u , 0, <^ ) and (/3 2 n , f3 l 3 n , B m , E m , s e , b). 
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Before we start to prove our main theorem, we need to make some preparations. 
Define W\ = s — sq, W 2 = /3 2 , W 3 = /3 3 , W4 = B — Bq, W$ = ip — ipo, we can derive the 
equations satisfied by Wj, i = 1, ■ ■ ■ , 5. 

(1 _ f^l)Q lWl + Q 2 W 2 + d 3 W 3 + \dx W 5 + diWi + d 4 W 4 + d 5 W 5 = -Fi(W, VW), 

<W 2 + + W 3 d 3 W 2 + ^Mf^d 2 W 1 - \d 2 W 5 + d 2 W 2 = F 2 (W, VW), 

u u 

dxW 3 + W 2 d 2 W 3 + W 3 d 3 W 3 + ^-^ds^! - \d 3 W 5 + d 3 W 3 = F 3 (W, VW), 

d{W 4 + W 2 d 2 W 4 + WjjSijWi = 0. 
AW 5 = p W! - eb(x) + F 5 (W, VW). 

Here 



(3.11) 



Fi(W, VW) = W 2 9 2 ^i + ^ 3 d 3 ^i " " ^r^Wi + ( A " "4 ^i^s 

- - a 2 W x + aiH^ + a!W 5 )d lSo + (jL - _ b 2 W 1 + 61 W 4 + hW 5 )d im , 

F 2 (W, VW) = C ^W 2dlWl - ^ dl W 5 - - C ^)d 2Wl 

Wq 

F 3 (W, VW) = C ^W 3 d lWl - ^ dl W 5 - - 



1 "1 "1 "0 



+ ( A - + - ^)^3 - (~ 2 ~ ^W^S, 



1 

F 5 {W,VW)=p-po-p W 1 . 



(3.12) 



and 



2c 2 (po) ( 7 -l) C 2 (po) 2c 4 ( P0 ) 
oi — 02 — 4 — , a 2 — g 1 1 — ) 

Un Ur, Un. 



~fc 2 (p )E 





61 = -j, di = -a 2 3is + & 2 <9i<A) = - 



c 2 (po) -^o) n cT 



1, c 2 (p ). , #0 



(3.13) 



d2 = d 3 = (^di^n o - ^i s o) 

d 4 = d 5 = ai9is - &i<9i<A) 



^0 (c 2 (po)-ng)' 

2E 



(c 2 (p ) -u 2 )u 2 

We emphasize that Fj((W, VW), i = 1, 2, 3, 5 are second order terms of Wi,i = 1, • • • ,5 
and does not contain any space derivative of W 2 , W 3 , W 4 . 
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Then it is easy to derive the elliptic system satisfied by W\ and W$: 



I Mil , -ikW-i - <h,Wr, ) + ( (^|^ - l)W 2 d 1 W 1 + ^-k^d 2 W l 



Un 



+d 3 ( _ i)w 3 aiWi + ^Wi) + d 2 f _ _ dlWl _ d5 w 5 



Un 



% + )« = Gi(W,VW), 











AW 5 = poW! - eb(x) + F 5 (W, VW). 
(^T^ " 1)^1^ - d 1 W 1 - cfe^ (0, x 2 , s 3 ) = *i(W, VW) +e(d2& + W + ^ . 



d 1 W 5 (0,x 2 ,x 3 ) = eE in (x 2 ,x 3 ), 
Wl(l,X 2 ,X 3 ) = es e (x 2l x 3 ), 
W 5 (l,x 2 ,x 3 ) = 0. 



(3.14) 



Here 



d(W, VW) =Si(Fi(W, VW) + d A W A ) + d 2 F 2 (W, VW) + d 3 F 3 (W, VW) 



+ d 2 ^iWiW 2 + d 5 ww 5 + ^2 w 2 ai^ 5 + w 2 Fi(w, vw) + d 4 w 2 w 4 j 

+ d 3 (diWiW 3 + d 5 W 3 W 5 + W 3 diVF 5 + W 3 Fi(W, VW) + d 5 W 2 W { 



+ 



2(d 2 W 3 d 3 W 2 - d 2 W 2 8 3 W 3 ) + d 2 d 4 W 4 + d2Fi(W, VW) 

1 



e6(x) + F 5 (W,VW) 



And W 2 ,W 3 and W4 will satisfy the following hyperbolic equations respectively: 



(3.15) 



d x w 2 + w 2 d 2 w 2 + w 3 d 3 w 2 + d 2 w 2 + 



( -k^d 2 W l - \d 2 W 5 = F 2 {W, VW), 



aiW 3 + w 2 d 2 ^ 3 + w 3 d 3 w 3 + d 3 w 3 + ^4^<Wi - -^3,w 5 = f 3 (w, vw), 

W 2 (0,x 2 ,x 3 ) = e(3 20 (x 2 ,x 3 ), 
W 3 (Q,x 2 ,x 3 ) = e/3 30 (x 2 ,x 3 ). 



(3.16) 



5iW 4 + W 2 d 2 W A + W 3 5 3 VF 4 = 0, 
W i (0,x 2 ,x 3 ) = eB in (x 2 ,x 3 ). 



(3.17) 
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3.3 Proof of Theorem I3TT1 



The main idea is simple: we construct an operator A on a suitable space, which will be 
bounded in a high order norm and contraction in a low order norm. 
The solution class will be given by 



W= (W 1 ,W 2 ,W 3 ,W 4 ,W 5 ) : ||W|| S = ||Wi,W £ 



5\\C a ^(Q e ) + ll^llc 2 '«(n e ) 

1=2 



< 5. 



Here 5 will be determined later. For a given W = (Wi, W2, W3, W4, W5) G H, we define an 
operator A : W = (Wi, W 2 , W 3 , W 4 , W s ) \-> W = (Wi, W 2 , W 3 , W 4 , W B ) mapping from S to 
itself, through the following iteration scheme. 
Step 1. To resolve W 4 . 



<9iW 4 + W 2 d 2 W 4 + w 3 d 3 W4 = 0, 
W 4 (0,x 2 ,x 3 ) = eB in (x 2 ,x 3 ). 



(3.18) 



The particle path (r, x 2 (r; x), :e 3 (t; x)) through (x±, x 2 , x 3 ), is defined by the following ordi- 
nary differential equations: 



cIx 2 (t;x.) 

dr 
dx 3 (r:x) 



dr 

x 2 (xi;x) = x 2 , 

5 2 (»i;x) = x 3 . 



W 2 (t, x 2 (r;x),5 3 (r;x)), 
W 3 (r, x 2 (t;x),x 3 (t;x)), 



(3.19) 



Since (W 2 ,W 3 ) G C^ Q (fi e ), it is easy to prove that (x 2 (0; x), x 3 (0; x)) belong to C AQ (ft e ) 
with respect to x. Indeed, the following estimate holds: 

||(x 2 (0;x),f 3 (0;x))|| c2 , Q(Ce) < C x . (3.20) 

One can find a unique solution W 4 (:z) = eB in (S: 2 (0; x), x 3 (0; x)) G C 1,a (ft e ) to ([53H]> . 
Furthermore, the following estimate holds: 



4||c 2 ><*(n e ) 



|C 2 > Q (T 2 ) 



< C 2 e. 



(3.21) 



Step 2. To obtain W\ and W5 by solving the following linearized second order elliptic 
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systems: 



<9i ( - 1)^1^1 - diWi - d 5 W 5 ) +d 2 [ - \)W 2 diW 1 + ^^-d 2 W l 

+d 3 ((^k^- - 1)# 3 £W + ^-d 3 W 1 ) + d 2 ((^^ - l)d 1 W 1 - d 1 Wi - cfe^s 
-^A)Wi- f^+5i(^)Vw5 = Gi(W,VW), 
AW 5 = PoV^i - eb(x) + F 5 (W, VW). 

u o / 
+e(d 2 /3? + d 3 p™ + \& 

d 1 W 5 (0,x 2 ,x 3 ) = eE in (x 2 ,x 3 ), 
W 1 (l,x 2 ,x 3 ) = es e (x 2 ,x 3 ), 
W 5 (l,x 2 ,x 3 ) = 0. 

(3.22) 

Take U\ = W\ — es e (x 2 , x 3 ), U5 = W5, then U\, U5 satisfy the following elliptic systems: 

di ((^k^ - 1)0^ - d 1 U 1 - d 5 U 5 ) + d 2 ( _ ^W^U, + ^%W X 

+d 3 ((^k^ - l)W 3 diU! + ^%Wi) + d 2 ( - 1)5! C/i - dil/i - 4C/5 

-^Po^i- f^+^(^)V^ 5 = Gi(W,VW), 
ti \« u y 

A<7 5 =ftf/i + G 5 (W,VW). 
( cJpo) _ l)diUi _ ^ _ d5Us \ 0jX2 ^ = h±(W, VW), 

d 1 U 5 (0,x 2 ,x 3 ) = H 2 (W,VW) = eE in (x 2 , x 3 ), 
Ui(l,x 2 ,x 3 ) = 0, 
U 5 (l,x 2 ,x 3 ) = 0. 



(3.23) 
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Here 



Gi(W, VW) = Gi(W, VW) + edxidxse) - V e^(-^^e) + e(cW 2 s e + ^U)- 
G 5 (W, VW) = -eb(x) + F 5 (W, VW) + e Po s e 



i=2 U U ° 



fTi(W, VW) = Fi(W, VW) + d 4 ^4 + e fyffi 1 + d^f + -^£ m + ed lS 



1 







(3.24) 

Set M = {U(x) G F^Oe) : {7(l,x 2 ,x 3 ) = 0.}. The weak formulation for (ET23]l is the 
following: To find U = (Ui, U5) eix M, such that for any V = (Vi, V 5 ) £ H x M, we have 
B(U,V) = #(V). Here 

B(U, V) = jT I ((^p " l)5if/i - dif/i - tWs) ^1 
+ ?(^-l ) ^ 1 C/ 1 + ^)^ 1 > 2 F 1 

V Uf, J 



(3.25) 



+ f (^M - 1)^^^ + ^^^L 

- ( da ((^^ - - dit/i - 4^5) V x dx 

+ / VU 5 -VV 5 dx+ / po^i^sdx. 



ff(V) = I GxV x dx- [ G 2 V 5 dx- I (H 1 V 1 + H 2 V 5 )(0,x 2 ,x 3 )dx 2 dx 3 . (3.26) 
Jn e Jn e Jt 2 

A direct computation shows that B is a bounded bilinear operator in EI x EI and $ is a 
bounded linear operator in EI x EI. Also B is coercive. Hence same arguments as in Chapter 
8 in shows us that Fredholm alternative theorem holds for ()3.23p . Hence uniqueness 
theorem will imply the existence of solution to (|3.23[) . It suffices to show that B(U, U) = 
implies U = 0. 

B(U,U) = / {(^-iXd^ + ^lid^f + idsUxf + lVU^Adx 

-[ [d 1 + d 2 (?-k^-l)}U 1 d 1 U 1 dx- [ dsdi^Usdx (327) 
+ I \% + d 1 {\))d 1 UsU l dx+ ! I ^ + d x d 2 \uldx 

Ja e \K n ) j ne ) 

+ jf (^po + d 2 d 5 ^U l U 5 dx = I + Y,I i = 0. 
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Now we estimate these terms respectively. 

h<^^(po + d 2 d b fU 2 dx + r/||£/ 5 ||£2 (f , e) 

~ ~kL ^o + ^5) 2 ^ + r?Cn e ||V^||| 2( n e) (3.28) 



h > min - max (^WU^ 2 

xg[o,i] xe[o,i] 

milljj.gfo !] po ||rr ||2 ^OMlrr li 2 

> 9 ^1 faro ) - max ^ 2 ,n v 

4 / 1 

^l/^^^max^^llVC/xlll^ 

+ max(^)[||f/illi 2(f , e) + lltfsll V) + livc/x|ii 2(ne) + ||vc/ 5 ||^ (ne) ; 

Hence we have 



= / + V/i> / \VUi\ 2 + \VU b \ 2 dx-C(b)( max (^Ull W- 



+ max (§)[||C/ 1 || 2 L2(Qe) + ||C/ 5 || 2 i2(f , e) + WVU^l^ 



xe[o,i] «q 

, llvrrr 112 A , ( min *e[0,l] Po C n e 2 \ n TT „2 

^ ^ / \max xe[0A] 2 xe[o,l] / 1 e> 



(3.29) 



(3.30) 



xe[o,i] ug " ^ 2 xe[o,: 



'o 



- ^ X |Wi ' 2 + |wsi2dx - ( ™] ( 4 )5||vt/ill ' 2(Qe) 

+ ™ (§)[ll^lli 2( Q e ) + ll^||! 2( n e) + l|VC/i|li 2(ne) + ||VC/ 5 ||i 2(ne) ]) 

>^(liv^illi 2(f , e) + livc/ 5 ||| 2(f , e) ). 

(3.31) 

Here we have used some special properties (12.220 of our background solutions. 

This implies that U± = U§ = 0. We establish the uniqueness for (|3.23p . By Fredholm 
alternative theorem, there exists a unique solution (Wi, W5) to (|3.22j) . 

By standard elliptic estimates in [lj, (|3.22p has a unique solution (Wi,Ws) G C 3,a (f2) x 
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C 3 ' Q (f2) and satisfies the following estimate: 
\\Wi,W 5 \\ c3 , a{ne) < CsAlGxCW, VW)|| c ^ (fie) + 6116(^)11^^) + ||^ 5 (W, VW)!!^^ 
+ 4d 2 (3T + d 3 f3 l 3 n + -^E in || C 2,a (T 2) + ||Fi(W, VW) + d4W4||oa.«(T») 

+ e||S m ||c2.«(T2) + e ll s ellc 3 ."(T 2 ) 

<C 4 (<5 2 + e). 

(3.32) 

Here C 3 depends only on background solution and O e , C4 depends also on /3™,f3™,B in and 
Step 3. To obtain W2, W-j by solving the following hyperbolic equations: 

8{W 2 + W 2 8 2 W 2 + w 3 d 3 w 2 + d 2 W 2 + ^-^d 2 ^i - \d 2 W 5 = F 2 (W, VW), 

d x W 3 + + W 3 d 3 W 3 + tfeW 3 + ~^-d 3 Wi - \d 3 W 5 = F 3 (W, VW), 

u o u o (3.33) 

W 2 (0, rE 2 ,x 3 ) = e/3 l 2 n (x 2 ,x 3 ), 
^ ^3(0,^x3) = e/3| n (x 2 ^3). 

Then by the characteristic methods, we have the following formulas: 
W 2 (x) = ee~ X 1 d2 ^ ds /3 l 2 n (x 2 (0; x), x 3 (0; x)) + jH e~ Z 1 d ^ ds (f 2 (W, VW) 

-^^8^ + \d 2 W 5 )(T,X 2 (T-^),X 3 (T^))dT, 

u o u o / (3.34) 

VF 3 (x) = ee-J'o 1 *W« i «^(x 2 (0;x),S3(0;x)) + £\~ ^ da ^ ds (f 3 (W ,VW) 

'"''^Wi + -i^WsVr.iaCr; x), x 3 (r; x))dr. 



(3.35) 



These enable one to obtain the following estimate: 

||(W 2 (x),W 3 (x))|| c2>ra(ne) < C 5 e||(/3r,/3f)|| c2 , ra(T2) + ||W||1 + W 5 ||c3.-(n e ) 

<C 6 (5 2 + e). 

This, together with the estimates (|3.32|) . (|3.35[) and (|3.2ip . gives 

||W|| a = ||^i,W 5 || C 3, a(ne) + \\W 2 ,W 3 , W4 C 2, a(ne) <C 7 (5 2 + e). (3.36) 

Here C7 = max{C 2 , C4, Ce}, which depends only on background solution and /3 2 n , /3 3 n , £> m 
and E in ,b,s e . 
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Choose ei small enough, such that if e < e±, then C 2 e < ^. Set 5 = 2Cje, then Cj5 < i 

and CV(<5 2 + e) < Cje +^5 = 5. This implies that A maps S to itself. 

It remains to show that the mapping A : E — > S is a contraction operator. Suppose 
A : (Wf,W$,Wg,w£,w!*) i ^ (W^W^W^W^Ws*) for fc = 1,2. Define the difference 
(Fx, y 2 , y 3 , n) = (^i 1 - VFi 2 , ^ - Wl Wi - Wl W\ - Wl Wl - Wg) and (Yl, Y 2 , Y 3 , Y A ) = 
{Wl - Wl W\ - wl wl - wl wl - wl Wl - Wg). 

Step 1. Estimate of I4. 

Indeed, Y4 satisfies the following equation: 



(3.37) 



3iY 4 + Wld 2 Y A + Wld 3 Y 4 = -Y 2 8 2 Wl - Y 3 d 3 Wi, 
Y 4 (0,x 2 ,x 3 ) = 0. 
Then the following estimate holds: 

< C 8 5||(y 2 ,y 3 )lbi.««(n e )- (3-38) 

Step 2. Estimates of Y\ and Y5. 

Y\ satisfies the following elliptic system: 

di ( - mr, - d lYl - d 5 Y 5 ) + d 2 ( (^?) _ + c ^ d2Yl ) 

+9 3 ( - IJWs 1 ^! + + d 2 ( _ 1)ftyi _ _ d5 Y 5 ) 

-k 0Y1 -(^ +dl{ ^) dlY5=HU 

AY 5 = Po Y 1 + F 5 (W\ VW 1 ) - F 5 (W 2 , VW 2 ), 

f(fM _ l)QiYi _ diYi _ d5 y\ (0)X2)X3 ) = Fl (W\ VW 1 ) - Fi(W 2 , VW 2 ), 
V u o J 

d 1 Y 5 (0,x 2 ,x 3 ) = 0, 

Yi(l,x 2 ,x 3 ) = 0, 

Y 5 (l,x 2 ,x 3 ) = 0. 

(3.39) 

h l = dew 1 , vwVdfw 2 , vw 2 )-a 2 ((^l-i)y 2 a 1 iy 1 2 'j _ 9s f(^°2 -i^w?) . 
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Then the following estimate holds: 



|(Yi,Y 5 )lb^ ( n e ) < C9(||ffi|| « ( n e ) + ||F 5 (W 1 ,VW 1 ) - F 5 (W 2 , VW 2 )|| cl , Q(Qe) 

+ ||Fi(w\ VW 1 ) -Fi(W 2 ,VW 2 ) +d 4 y 4 || c i,« (T 2 ) ) (3.40) 



< Cio^^||Yi,i > 5|| c -2, a ( ne ) + ||Y 2 , Y 3 , Y 4 || c i, Q (Q e )^. 



Step 3. Estimate of I2, Y3. 

It follows from (I3,33|) . Y 2 and Y 3 satisfy the following system: 



d 1 Y 2 + wid 2 Y 2 + w£d 3 Y 2 + d 2 Y 2 + ^%Wi - \d 2 Y 5 = k x , 

8iY 3 + WfaYa + W£d 3 Y 3 + d 3 Y 3 + - \d 3 Y 5 = K 2 . 



Y 2 {Q,x 2 ,x 3 ) = Q, 
Y 3 (0,x 2 ,x 3 ) = 0. 



Here 



K x = -(J 2 d 2 Wl + %d 3 Wl) + F2CW 1 , VW 1 ) - F 2 (W 2 , VW 2 ), 

K 2 = -(Y 2 d 2 Wi + %d 3 Wi) + F 3 (W 1 , VW 1 ) - F 3 (W 2 , VW 2 ). 
Then the following estimate holds: 



||(F 2 ,y 3 )||c^(n e ) ^C'nl <J||(y 2 ,y3,n)|| t 7i.a(n e ) + «l|Vyi,Vy 5 ||c7i.«(n e ) 

< c 12 s( ||yi,f'5|| f7 2,a( ne ) + ||y 2 ,y 3 ,y4|| c .i,«(n e ) ). 



(3.41) 



(3.42) 



(3.43) 



Setting C\ 3 = max{Cs, C\q, C\ 2 }, then take e 2 small enough such that Ci 3 e 2 < 1. Now 
we take eo = min{ei,e 2 }, then if < e < eo, A maps E to itself a nd is contraction in low 
order norm. Hence A has a unique fixed point, which will be the solution to (jl.ip . We have 
finished our proof. 

Remark 3.2. We can not generalize this result to a general 3-D nozzle. In general, the velocity 
field f3 2 and j3 3 will lose one order derivative when integrating along the particle path. The 
delicate nonlinear coupling between the hyperbolic modes and elliptic modes (/3 2 ,/3 3 ) makes it 
extremely difficult to develop an effective iteration scheme in a general 3-D nozzle. How to 
effectively explore the good property of the quantity W = <9 2 /3 3 — <9 3 /3 2 + /3 3 <9i/3 2 — P 2 d\f3 3 will 
be investigated in the forthcoming paper. 
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4 Appendix 



In this appendix, we give the detailed calculations for (|2.6p . First we rewrite (|2.22p as 
follows: 

G _1 D/3 2 + n/n ; 1 - - (feft - a 2 )(v9 - %)) = 0, 

(4.1) 



2(B + ^-%)) 
G _1 D / S3 + nrn , 1 ^^ (Mi - da)(<p - h{ P )) = 0. 



2(B + <p-h(p)) 

Applying f3 3 di — d 3 and —(j3 2 d\ — d 2 ) to (|4.ip and adding them together, we obtain 



= (Mi - 5 3 )(G- 1 D/3 2 ) - (Ml - 9 2 )(G- 1 D/3 3 ) 



+ 



2(B + <p-h(p)) 



(Ml - d 3 )(Mi ~ d 2 ) - (Mi - d 2 )(Mi - d 3 ) 
(Mi - d 3 )B(Mi ~ d 2 )(<p - h{p)) 



(<p - h(p)) 



2(B + ^p-h(p)f 

- (Mi - d 2 )B(Mi ~ 3 )(<p - h(p)) 

(Mi - 5 3 )(G- 1 D/3 2 ) - (Mi ~ 3 2 )(G- 1 D/3 3 ) H- 
1 



1 



2(B + if-h( P )y 



2(B + <p-h(p)) 
(l,ft,ft)-(V(^-Mp))xVB 



d 1 (<p-h(p))W 



While 



J : = (Ml ~ 3 3 )(G- 1 D/3 2 ) - (Mi ~ a 2 )(G" 1 D/3 3 ) 
= G- 1 D((/3 3 a 1 - a 3 )/3 2 - (Mi - 5 2 )/3 3 ) 

3 

+ G- 1 J2((Mi - Wjdjfo - (Mi - ddPjdjfo) 

3=1 

3 r 

+ G- 1 J2Pj\ ii^di - d 3 )d j - dj(Mi - «9 3 )]/3 2 
3=1 L 

- [(Mi - d 2 )dj - dj(Mi - o 2 )]f3 3 

+ (Mi - 3 3 )G- 1 D/3 2 - (Mi - 3 2 )G~ 1 D/3 3 
=G- 1 VW + Ji + J 2 + J 3 . 
Now we compute Jj, i = 1,2,3 respectively. 



Ji = G- 1 



w(d 2 /3 2 + d 3 f3 3 ) - (Mi - 3 3 )/W3 + (Mi - d 3 )M^2 



(Mi - d 2 )M2f3 3 + (Mi - d 2 )M2h 



G- l [w(d 2 p 2 + d 3 p 3 ) + z}. 



(4.2) 



(4.3) 



(4.4) 
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Here Z = fad x fad 3 fa - fadifad 3 fa + fad x fad 2 fa - fad!fad 2 fa. 

3 r 

■h = G- 1 Y,Pj\ - d 3 )d 3 - djifod! - d 3 )\fa 



hdi - d 2 )dj - djihdi - 8 2 m 



G~ l Z. 



Here 



Jz = -G 



-2 



(fab - d 3 )GT>fa - (fob - d 2 )GBfa 
3 



/?i(/3 3 3i - d 3 )faT>fa - frifodi ~ d 2 )pjBp 3 



= -2G~ 2 [W(faVfa + faDfa) + J31]. 



J31 = fa(fadi - d 2 ) faDfa + fa(fad! ~ d 3 )faBfa 

- P 3 (Mi ~ d 3 )faT>fa - W 2 di - d 2 )faT>fa 



(fed! - d 2 ) faDfa - {fab - d 2 ) faDfa 



+ /3 3 



(Mi - d 3 ) faDfa - (fad! - d 2 ) faDfa 



(4.5) 



(4.6) 



(4.7) 



= fa (fad! - d 2 )fad!fa - (fad! - d 2 )fad!fa 
+ fa {fad! - d 3 )fad x fa - (fad! - d 3 )fad!fa 
+ (fad! - d 2 )fa(pld 2 fa + fafad 3 fa) - (fad! - d 2 )fa(fcd 2 fa + fafad 3 fa) 
+ (fa^ - d 3 )fa(fafad 2 fa + (3ld 3 fa) - (fad! - d 3 )fa(fafad 2 fa + pld 3 fa) 

= GZ. 

Hence we have J 3 = -2G~ 2 [W(faDfa + faDfa)] - 2G~ 1 Z. Substitute these calculations 
into the above formula to get: 
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= J + 



-M<p-h{p))W 



2(B + <p-h(p)) 



2{B + ip-h{p)Y 



(l,/32,/9 3 ) • (V(<p-h{p)) xVB 



1 



1 



2(B + ^-/ l (p)) 2 
G _1 DW - G -1 WDs - G~ 2 VFDG 
1 



2(5 + y>-%)) 
(l,A,ft)'(v(^-%))xV5 

1 



P G' 2(B + <p-h(p)y 
This implies that 



2(B + <p -h{p)Y 
(l,fc,fc)-(v{ip-h{p))xVB) 



1 



pG 7 2p( J B + v ?-/i( / o)) 2 



(1,02,03)- fa(<p-h(j>)) x VS) 



0. 



(4.8) 
(4.9) 
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